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ABSTRACT 

We present Herschel-SPIKE Fourier Transform Spectrometer (FTS) and radio follow- 
up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed dis- 
tant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) 
we detect [Oiii]88/im and [Cii]158/xm lines at a signal-to-noise ratio of ~ 5. We 
do not have any positive line identification in the other fainter target H-ATLAS 
J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with 
positive line detections with the FTS, with continuum flux just below 200 mJy in the 
200-600 /im wavelength range. The derived redshift of SDP.81 from the two detections 
is z = 3.043 ± 0.012, in agreement with ground-based CO measurements. This is the 
first detection by Herschel of the [O ill] 88 /im line in a galaxy at redshift higher than 
0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation 
region (PDR) models with different physical conditions, we derive the PDR cloud 
density n w 2000 cm"'^ and the far-UV ionizing radiation field Go « 200 (in units 
of the Habing field - the local Galactic interstellar radiation field of 1.6 x 10~^ W 
m~^). Using the CO derived molecular mass and the PDR properties we estimate 
the effective radius of the emitting region to be 500-700 pc. These characteristics are 
typical for star-forming, high redshift galaxies. The radio observations indicate that 
SDP.81 deviates significantly from the local FIR/radio correlation, which hints that 
some fraction of the radio emission is coming from an AGN. The constraints on the 
source size from millimiter-wave observations put a very conservative upper limit of 
the possible AGN contribution to less than 33%. These indications, together with 
the high [O ill]/FIR ratio and the upper limit of [O i] 63 //m/[C il] 158 /xm suggest that 
some fraction of the ionizing radiation is likely to originate from an AGN. 

Key words: galaxies: evolution — infrared: galaxies — submillimetre: ISM — radio 
continuum: galaxies — galaxies: individual (SDP.81: H-ATLAS J090311.6-H003906, 
SDP.130: H-ATLAS J091305.0-005343) 



1 INTRODUCTION 

The detection of large numbers of high-redshift, optically 
obscured, massive star-forming galajcies by ground-based 
sub-mm instruments such as Submmilimetre Common-User 
Bolometer Array (SCUBA) and Max Planck Millimetre 
Bolomet er (MAMBO) was a surpris e and a major discov- 
ery (e.g. ISmail. Ivison fc Bla"inlll997^ . These submm galax- 
ies (SMGs) present amongst the most difhcult obstacles for 
models of galaxy formation and evolution (e.g. iBaugh et al.l 
l2005l ). The SMGs are a subset of the dusty star-forming 
galaxies that account for more than half of the far-infrared 
background (|Devlin et al.|[2009l '). Their star- formation rates 
are orders of magnitudes higher than most galaxies in the 
local Universe and they are considered the precu rsors of to- 
day's most massive elliptic al or bulge galaxies (|Lillv et al.l 
ll999l : ISwinbank et al.|[200^ ). 

Due to the fact that more than half of the light from the 
SMGs is absorbed and reprocessed by dust, it is exceedingly 
important to study these sources at their peak emission, 
which occurs in the far i nfrared (FIR). Prior to the Her- 
schel (|Pilbratt et al.|[201(]| '). ground based and balloon-borne 
studies of SMGs were limited to small fields and produced 
modest samples detected at low significance. This situation 
is changing dramatically with the number of wide-field sur- 
veys with Herschel Space Observatory, such as the Herschel 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
t E-mail:ivaltchanov@sciops.esa.int 



Astrophysical Terahertz Large Area Survey (H-ATLAS). H- 
ATLAS is the larg est Herschel Open Time Key Programme 
(|Eales et all l2010f l. and wiU map 550 deg in six extra- 
galactic fields in five photometric bands, u sing in parallel 
mode both the PAC S (jPoglitsch et al.ll201(i ) and the SPIRE 
(|GrifHn et al.ll2010l ) instruments. One of the most interest- 
ing results of H-ATLAS will be the detection of a signifi- 
cant number of gravitationally lensed objects in the distant 
{z > 2) Universe. This is possible in such a shallow but large 
area survey because of the combined effects of a strong neg- 
ative if-correction and steep number counts in the SPIRE 
500- ^m waveband (e.g. [Blain 199^). Objects brighter than 
100 mJy at 500 /im will be a mixture of lensed high-redshift 
galaxies, nearby galaxies and flat-spectrum radio sources 
(|De Zotti et~alll2005l : iNegrello et"alll2007l : lEales et al.ll2010l ) 
and current models predict that H-ATLAS will detect some 
^ 350 strongly lense d objects at 500 pm at redshifts > 1 
(|Negrello et al.ll2007l ). 

High-redshift, bright SMGs open a new possibility for 
the detailed study of the physical conditions of the inter- 
stellar medium, at rest-frame far-infrared. Atomic emis- 
sion lines in this range provide the major coolants for 
neutral gas exposed to the far-ultraviolet (FUV) ioniz- 
ing flux of nearby early-type stars. Star-forming molecu- 
lar clouds where the physical processes are dominated by 
FU V photons are called photo- dissociation regions (PDRs, 
e.g. iHollenbach fc Tielenslll997l ). FIR lines are not strongly 
affected by extinction in these regions, and so provide an 
excellent probe of the embedded stellar radiation fields 
and their effects on the physical conditions of the neutral 
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gas. The majority of their cooling o ccurs via atomic fine 
structure hues of [Cll], [Ol], [Olll] jBrauher et all 120081') 



and the molecular rotational ^'^CO lines ( Kramer et alj 



I2OO4I '). There have been several pie- Herschel efforts to study 



galaxies, such aslStacev et alJ (Il99ll. 


2010l):lMalhotra et all 


1I2OOII): Maiolino et al. (|2005|, 2005 


); 


Brauher et al. (l2008h: 


1 Walter et al. (l2009a"b');'WaEE et al 


(I2OIGI): Ferkinhoff et all 


l2010l):IVasta ct al. (20f0); iHailev-Dunsheath et al.l (I2OIOI). 



as wel l as molecular (CO) studies o f SMGs, e.g. Greve et al.l 



ll2005h:lTacconi et al.l ||2006|. l2008h: iDaddi et all (|2009al lbt): 
iRiechers et all (|201G| '): iDanielson et all (120111 ). Herschel, 
however opens up a new frontier with its highly sensitive 
spectrometers, especially for observations of lines like [C 11] 
at 158 fim at redshifts up to ~ 3. 

During Herschel's Science Demonstration Phase, one 
of the H-ATLAS fields, GAMA9, w as o bserved, covering 
~ 4 X 4 deg^ (see lSmith et all 1201 ll and lRigbv et al.ll201ll 
for the H-ATLAS data in this field) . The proce s sing of the 
Herschel observati ons is presented in Ibar et ahl (|2010bl ) for 
the PACS and in iPascale et al] (|201ll ) for the SPIRE ob- 
servations. We have identified 11 sources with flux density 
at 500 fim S500 ^ 100 mjy and the results of the analy- 
sis of these candidate lensed sources, from follow-up obser- 
vations using facilities both on the g round and in space , 
are presente d in Negrello et al. J2OIO') ; 'L upu et"al] ([2O1O) ; 
iFraver et all (|2011^ : iHopwood et al, (.201 llV 

The subject of this paper is to present the SPIRE 
Fourier- Transform Spectrometer (FTS) and radio obser- 
vations and report on the first r esults from the far-IR 
spectroscopy for two targets from iNegrello et al.l (|201(]| ): 
H-ATLAS J090311.6-h003906 designated SDP.81 and H- 
ATLAS J091305. 0-005343 designated SDP.130. Their po- 
sitions, measured redshift and photometric properties in 
the FIR are shown in Table [1] Each of these two tar- 
gets was initially associated with an elliptical galaxy, yet 
based on their Spectral Energy Distribution (SED), peak- 
ing at ~ 350 i-im, it was clear that the foreground galaxy 
acts as a lens for a background source. Subsequent follow- 
up using Caltech Submmilimeter Observatory/Z-Spec and 
Green Bank Telescope/Zpectrometer instruments detected 
several CO lines at high signal-to-noise ratio (SNR) and se- 
cured a spectroscopic redshift of 3.04 2 ± 0.001 for SDP.81 
and 2.625 ± 0.001 for SDP.130 (s ee INegrello et al] I2OI0I : 
ILupu et al]|2010l : iFraver et al]|201ll for more details). 

The feasibility of the proposed FTS follow-up was 
well demonstrated during SPIRE's Performance and Ver- 
ification phase. One bright, lens ed, high-redshift SMG, 
SMM J2135-0102 at 2 = 2.3 (jSwinbank et al.l |2010| ). 
was observed with the SPIRE FTS, achieving a de- 
tection of the [C 11] 158/xm line at 4.3(t (line flux of 
(1.7 ±0.4) X 10"^'' W m ■^) as well as tentati ve detections 
of [O l] at 145.5/im and [Nil] at 122.1/im (|lvison et al.l 
l2010dl ). 

The structure of the paper is as follows: in Section [2] 
we present the FTS and radio observations and the FTS 
observing strategy. Some of the characteristics of the en- 
vironment (telescope and instrument) during the FTS ob- 
servations, relevant for the subsequent data processing are 
also presented. The FTS data processing, with some details 
on the deviations from the standard Herschel-SPIKE FTS 
data processing, is described in Section (3) The results are 



presented in Section |4] and summarised in the conclusions 
(Section [5]). 

For all cosmological dependent parameters we use h — 
0.73, »m = 0.24 and fl at geometry, i.e. f^A = 1 - n„ = 0.76 
(|Lahav fc Liddlell2006l ). With these parameters the luminos- 
ity distances to redshifts 2.625 and 3.042 are 22.06 and 26.41 
Gpc, respectively. 



2 OBSERVATIONS 
2.1 SPIRE FTS 

The FTS observations were carried out on June 1st 2010 
{Herschel observing day number 383). We targeted the two 
H-ATLAS lensed candidates, SDP.81 and SDP.130 and a 
nearby (offset at ~ 2') dark sky region, devoid of sources and 
chosen using the available SPIRE maps. The three observa- 
tions were done in a sequence: SDP.81 then the dark sky 
and then SDP.130, each one with exactly the same SPIRE- 
FTS observing template: 100 repetitions using single point- 
ing mode, sparse spatial sampling and high spectral reso- 
lution. The total time for each observation, including over- 
heads, was 3*^52.2™, with 66.6 s on-source time per FTS scan 
(one repetition consists of one forward and one reverse scan, 
i.e. two scans per repetition). 

T he measured absolut e pointing error of Herschel is 2" 
(68%, IPilbratt et ai]l2010l ) - this is much smaller than the 
smallest FTS beam of ~ 17" at 200 nm (|Swinvard et all 
|201(]| ). The pointing stability of Herschel during the obser- 
vation was excellent: during the 4 hours of spectral scans 
the telescope jitter stayed within 0.14" (68%) from the tar- 
get position. The largest deviation was 0.54". 

The main source of emission in the FTS bands is the 
Herschel telescope itself, although this is a constant back- 
ground as the telescope temperature does not change signif- 
icantly on the time scale of our observations. The emission 
of the telescope ranges from 100 Jy to 1000 Jy in the FTS 
bands (Swinyard et al., in prepration). Targets like SDP.81 
and SDP.130, with continuum levels around or below 200 
mJy are just a tiny fraction of this signal. That is why it 
is very important to have an off-source reference dark sky 
measurement, performed as close in time as possible to the 
target observation, in order to be used to subtract the tele- 
scope and instrument contribution from the source spectrum 
(see Section |3]). 

By design, the FTS also has a second input port that 
sees a spectrometer calibration source (SCAL) and this is 
the second source of significant thermal emission. Regard- 
less of t he fact that SCAL is not used and the lamp is turned 
off (see lSPIRE OMll201Cl^h . it stays at the instrument enclo- 
sure temperature and emits thermally in the two FTS bands. 
This instrument contribution is also removed with the help 
of the dark sky observation. As SDP.130 is about 30% fainter 
than SDP.81 at 350 ^im and at a different redshift (see Ta- 
ble [TJ, we have added this observation as a second dark sky 
measurement in order to decrease the noise of the reference 
subtraction for SDP.81. 



^ The SPIRE Observers' Manual is available from the Herschel 
Science Centre: 

|http : //herschel ■ esac . esa. Int/Docs/SPIRE/pdf /splre_om.pdf | 



4 Valtchanov et al. 



Table 1. Positions and photometry measurements of SDP.81 and SDP. 130. The quoted positions are those determined using IRAM 
PdBI at 0. 62" X 0.30" beam reso lution (Neri et al. in preparation), thus accurate to < l". The redshifts and the SPIRE photometry are 
taken from iNegrello et al. while the PACS 160 /im measurements are from this work. The upper limits at 70 /^m arc Icr. 



Target ID 


Rcdshift 


RA 


Dec. 


■570 


Sim 


■S'250 


■Saso 


S500 






{J2000) 


{J2000) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


SDP.81 


3.042 ±0.001 


09:03:11.6 


+00:39:06.7 


< 8 


51.4 ±5.0 


129.0 ±6.6 


181.9 ± 7.3 


165.9 ±9.3 


SDP. 130 


2.625 ±0.001 


09:13:05.3 


-00:53:42.8 


< 9 


45.3 ±7.7 


105.0 ±6.5 


127.6 ± 7.2 


107.5 ±9.0 



During the observations the primary mirror tempera- 
ture was on average 88.20 ± 0.06 K, the secondary mirror 
temperature stayed at 84.290 ± 0.004 K with variations at 
the mK level. The instrument enclosure temperature dur- 
ing the SDP.81 observation was 4.70-4.75 K, during the ofT- 
source dark sky measurement it was 15 mK colder and 22 
mK colder during the SDP. 130 observation. Such seemingly 
small differences in the instrument and telescope tempera- 
tures, even at a level of a few mK, during the dark sky and 
the target measurements, can lead to overall additive e ffect 
on the continuum level. Following llvison et al.l ([ioiO^, we 
have used the three SPIRE photometry points at 250, 350 
and 500 /im in order to recover the continuum level (see 
Section [S]) . 

2.2 PACS photometry 

In addition to the SPIRE FTS spectral scans we did short 
PACS cross-linked mini-scan maps for each source at 70 and 
160 ^m (PACS blue and red bands). This was necessary be- 
cause SDP.81 was not bright enough to be detected in the 
PACS parallel mode scans, while SDP. 130 was not in the re- 
gion covered by PACS in the PACS/SPIRE parallel mode. 
The PACS mini-scan maps were processed with version 5.0 
of the Herschel Interactive Processing Environment (HIPE, 
lOtt et al.|[2O10l ) using a standard data reduction script ap- 
plying a temporal high-pass filter to remove the 1/f noise 
after masking the source for the filtering. The data cube was 
projected with the photProject task and the 2 scan directions 
were combined with the Mosaic task. Unfortunately the cho- 
sen integration time of ~ 15 minutes was not sufficient to 
detect the two sources at 70 fim, where the uncertainties in 
the SED were the highest. The PACS flux density of both 
targets was measured in an aperture of 10" radius, corrected 
for the encircled energy fraction (see the PACS Observers' 
Manualfl) and are reported in Table [T] At 70 /im we only 
report the la upper limits. 

2.3 Radio observations 

SDP.81 was observed by the Extended Very Large Ar- 
ray (EVLA) for two 3-hr tracks during 2010 July 17-18 
(programme AI143/10A-253), when the EVLA was in its 
most compact configuration (D), recording data every Is 
and using 16 Ka-band antennas. Antenna pointing was 
checked at 5 GHz. Much of the first track suffered from 

^ The PACS Observers' Manual is available from the Herschel 
Science Centre: 

|http : //herschel ■ esac . esa. Int/Docs/PACS/pdf /pacs_om.pdf] 



a correlator problem and was discarded. For SDP.81, the 
^^CO J = 1-0 line (CO(l-O) hereafter) is redshifted to 
28.53 GHz for z = 3.042, a frequency at which just one 
of the two sub-band pairs of the EVLA's new Ka-band re- 
ceivers could be utilised, yielding 128 MHz of instantaneous, 
dual-polarisation bandwidth during the earliest shared-risk 
phase of EVLA commissioning with the new Wide band In- 
terferometric Digital Architecture (WIDAR) correlator. For 
our observations at 8.40 GHz, we obtained 256 MHz of con- 
tiguous, dual-polarisation bandwidth. In total, we recorded 
2 hr of data at 8.40 GHz and a similar amount at 28.53 GHz, 
resulting in an overall on-source integration time of ~ 85 
minutes. 3C 286 was used for absolute fiux calibration, and 
the local calibrator, J0909+0121, was observed every 5min. 

Editing and calibration were accomplished within 
AIPS (31DEC10). For 3C286, we used an appropri- 
ately scaled 22.5-GHz model to determine gain solutions; 
J0909+0121 was used to determine the spectral variation 
of the gain solutions (the "bandpass"), after first remov- 
ing atmospheric phase drifts on a timescale of 12 s via self- 
calibration with a simple point-source model. We were able 
to employ standard AIVS recipes throughout the data- 
reduction p r ocess, following the techniques employed by 
llvison et al.l (|2010bl ). The synthesized beam for the CO ob- 
servations is 3.3" X 2.4" and « 12" x 8" for the continuum. 

Figure [T] summarizes the result of the EVLA observa- 
tions for SDP.81. For CO(l-O), we measure an integrated line 
flux of 1.26 ±0.20 Jykm s~^, consistent wit h the GET value 
of 1.11 ±0.25 Jykms-i (jPraver et al.ll201lD . the Full Width 
at Zero Intensity (FWZI) is ~ 750 kms"\ At 8.40 GHz, we 
measure a total flux density of >S'8.4GHz = 700 ± 50 /xJy. A 
contour plot of the integrated CO (1-0) emission is shown in 
the inset of Figure [T] overlaid on a greyscale 880 fim contin- 
uum emission from the Submillimeter Array (SMA). 

SDP.81 was observed and detected with the Institut 
de Radioastronomie Millimetrique (IRAM) Plateau de Bure 
Interferometer (PdBI) at 240 GHz and at resolution of 
0.62" X 0.30" on 2010 February 15th (Neri et al. in prepara- 
tion). SDP.81 was also observed at higher resolution with the 
Combined Array for Research in Millimeter-wave Astron- 
omy (CARMA) on 2010 February 22nd and 23rd (program 
cx306), when CARMA was in its most extended A config- 
uration at 240 GHz. Due to previously arranged scheduling 
plans, the phase monitoring array was not deployed dur- 
ing these observations. However, the phase rms was less 
than 100 /im and the atmospheric opacity at 230 GHz (or 
T"230GHz) was less than 0.2 at zenith, allowing the observa- 
tions to go ahead without the monitoring array. The FWHM 
size of the test calibrator was 0.19" x 0.14", while the theo- 
retical beam size is 0.17" x 0.14". This is pretty good agree- 
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Figure 1. Spectrum of ^^CO J = 1-0 line emission from 
SDP.81 from the Expanded Very Large Array. Inset: Contours of 
tlie CO(l-O) emission, integrated across the velocity range —400 
to +350 km s"'^ (indicated by the line below the spectrum, and 
where v = Okms"! at 28.51836 GHz for z = 3.042), which corre- 
sponds roughly to the FWZI of the line, overlaid upon a grayscale 
S MA 880 lira con t inuum image (scaled from 0-8 mjy) presented 
by iNegrello et al ] l|20ld '). The synthesized beam for the CO ob- 
servations is shown in the inset, bottom left (3.3" X 2.4"), and 
the contours are plotted at —2.5, 2.5,3.5... 6.5 Xc, where the rms 
noise level, a = 100 /ijy beam" ^ . The synthesized beam for the 
continuum observations is ?s 12" X 8" and the size of the image 
is 20" X 20". 

ment and suggests atmospheric "seeing" was not a big prob- 
lem during the observations. 

The CARMA observations of SDP.81 were in two ~ 5- 
hr tracks, the noise obtained from the first track being 
much higher than the second, due to the higher atmospheric 
opacity. Therefore, the two tracks were not coadded and 
only the results from second track that was observed in 
T"230GHz = 0.1 are presented here. The data were reduced 
with and without a gaussian taper of 0.3". The taper reduces 
the resolution (and sensitivity), but may pick out low- level 
emission. SDP.81 was not detected in either case, and the 
240 GHz respective resolutions and sensitivities reached with 
CARMA were 0.17" x 0.14" and 0.67mJy/beam in the no- 
taper and 0.30" x 0.26" and 1.1 mjy/beam in the 0.3"-taper 
maps. 

SDP.81 is also present in the FIRST survey 

jBecker et all 1 19951 ) at 1.4 GHz with 5'i.4GHz = 950 ± 
150 /im. 



3 SPIRE FTS DATA PROCESSING 

In this section we briefly describe our data processing ap- 
proach and underline where it deviates from the standard 
SPIRE FTS data processing. Detailed description of the 
standard FTS pipeline can be consulted in the SPIRE Data 
Reduction Guidelf] 

^ Available at I http: //herschel ■ esac ■ esa. Int | 



Each of the three observations, two targets and one ref- 
erence dark/off-target, are processed separately out to the 
final spectrum. The forward and the reverse scans are kept 
separated throughout the whole processing. We consider 
only the central detectors of the two SPIRE FTS short- 
wavelength (SSW) and long-wavelength (SLW) bolometer 
arrays. Note that the off-axis bolometers are of a limited 
use as the optical path and the standing-wave (fringe) pat- 
tern in the interferograms are very different. The fringes are 
at high optical path differences and consequently they have 
the undesirable effect of adding a fixed pattern noise that 
affects the high resolution part of the spectrum. 

The FTS detector timelines, sampled at 80 Hz, are first 
cleaned of glitches due to high energy particles hitting the 
substrate, the bolometers or the readout electronics. The 
glitch removal utilises a wavelet based method. The fraction 
of samples identified as glitches was ~ 0.5% for the SSW 
and ~ 0.8% for the SLW for the three observations. Next, 
instrumental effects are removed: electrical cross-talk and 
time domain phase correction. Then we convert the bolome- 
ter recorded voltage into its temperature using a calibra- 
tion between the bolometer resistance and temperature, es- 
tablished during ground testing. This is in contrast to the 
standard voltage-based pipeline FTS processing, where a 
non-linearity correction is applied to the read-out voltages. 
The reason for using bolometer net temperature is that it 
represents much more closely the physical response of the 
bolometer to incident radiation and therefore "naturally" 
accounts for any non-linearity and drift in that response. In 
the circumstances we are faced with here, where we are crit- 
ically dependent on being able to accurately subtract two 
large signals to give a very small one, we have found that 
this, combined with the interactive phase correction steps, 
gives a small but significant improvement over the current 
pipeline processing. It may be that using the temperature 
based pipeline ends up being unnecessary with further im- 
provements in the voltage based pipeline, but for the time 
being we prefer this method, to remove any dependency on 
the drift correction and non-linearity steps which could in- 
ject systematic uncertainties into the data. 

The resulting detector timelines are then converted 
to interferograms with the use of the spectrometer mirror 
mechanism (SMEC) timeline that provides the optical path 
difference (OPD) at each time sample. Thus, at the end of 
this step we have interferograms of signal versus OPD. The 
100 forward and 100 reverse interferograms are then divided 
into 2^* parts, which we shift iteratively in zero-path differ- 
ence (ZPD) until the interferograms are precisely aligned, 
removing any residual phase shift introduced by the detec- 
tors' temporal response and beamsplitters jSpencer et al.l 
l20ld : iNavlor et all [2OI0I ') . The interferograms from all 200 
scans are then merged, the forward and the reverse ones 
separately, they are baseline corrected and during the aver- 
aging process the data points which are marked as glitches 
are discarded. The next steps include gain correction and 
phase correction of the interferograms, which are further 
zero padded and converted to the spectral domain using Fast 
Fourier Transforms. To convert the flux to physical units 
(Jy) and to correct for the relative spectral response func- 
tion (RSRF) we use a deep observation of Uranus. At this 
stage, the forward and reverse interferograms are averaged 
together. This procedure was followed for the two targets 
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Observed frequency (GHz) Observed frequency (GHz) 

Figure 2. The FTS high and low resolution spectra of SDP.Sf (left) and SDP.fSO (right) overlaid with the three SPIRE photometric 
points (dots) and the SED template fit. The photometric errors, even including the 7% flux calibration uncertainty of the SPIRE 
Photometer, are smaller than the dots. The expected positions of a number of fine-structure lines are indicated with vertical bars. 



and the reference sky. The last step is the subtraction of the 
reference sky from the target and the resulting spectra of 
SDP.SI and SDP.130 are shown in Figured 

As we described in the previous section, the absolute 
continuum calibration, especially for faint targets, is quite 
difficult. Small changes, of the order of few tens of mK, in 
the temperature of the instrument and the telescope during 
the reference dark sky observation and the target observa- 
tions can lead to quite significant offsets in the continuum 
level. We iteratively adjusted the overall dark sky level in 
order to get the best match with the available photometric 
points (Swinyard et al., in preparation). A further correction 
was applied to the SLW spectra for both targets in order to 
remove the effect of variations in instrument temperature 
ISpencer et al.ll20f(]| ). These corrections are already applied 
to the spectra shown in Figure [21 Note that the continuum 
offset affects the low resolution continuum too. This makes 
the use of the FTS in low resolution mode, for an SED-like 
measurements, not very suitable, especially for faint targets. 



however we did not perform this additional smoothing as 
we are looking for very faint features. The SPIRE FTS in- 
strumental line shape is very close to a Sine function, with 
noticeable asymmetry of the first high-frequency minimum. 
The reason for this asymmetry is still under investigation 
by the SPIRE instrument team. This may have very small 
effect on the line flux measurements, however the line cen- 
troids are not affected, as it was shown that the centroids 
of CO ladder lines in n earby sources are as good as 1/20 of 
the resolution element (|SPIRE OMII2OIOI ). 

Thanks to the already available h igh accuracy spectro- 
scopic redshifts f o r the two galaxies iNegrello et al.l 120101 : 
ILupu et all I2OIOI : iFraver et all I2OIII ). we know where to 
search for lines at specific wavelengths. We fit a Sine 
function model to the continuum subtracted spectrum in 
±5000 kms~^ rest-frame, around the expected positions of 
a number of ISM lines: 



I{v) = /(O) 



sin{{v — vo)/va) 



(1) 



4 RESULTS 

The FTS covers the wavelength region from 197 to 670 /im. 
Some of the most prominent and well studied ISM cooling 
lines fall in this range, for the redshifts of the two targets 
in this study. The FTS spectral resolution is l^v = 1.2 GHz 
for the high resolution mode and it is constant in frequency 
space. This spectral resolution corresponds to instrumental 
line FWHM from 280 to 970 km s"\ over the FTS spectral 
range. The measured FWH M of the CO(l-O) line is 440 ±50 
kms~^ iFraver et al.ll201ll ). also compatible with the EVLA 
measurement (see Figure [!}. This means that some lines 
could be marginally resolved at A < 300 ^m. 

The instrumental line shape of all FTS instruments is 
a Sine function due to the truncation of the interferogram 
by the limited travel of the moving mirror. This function 
is characterised by one central peak and a number of sec- 
ondary maxima at lower intensity. There are techniques to 
suppress the secondary maxima at t he price of lowering th e 
resolution, the so called apodization (|Navlor fc Tahicll2007t ). 



where Va is the Sine line width in kms^^, measured be- 
tween the two zero crossings nearest to the peak and no 
is the line centroid offset from the expected position for 
the given redshift. For the high resolution mode we have 
Vcr = 1.193 (Ac//im) kms~ at the line centroid Ac, we keep 
this parameter fixed. Note that the line FWHM is ufwhm = 
Wct/1. 20671. The output parameters are the line peak 1(0) 
and the line centre offset wo- The integrated line flux for for 
a Sine function is then calculated using: 

5'iinc = 0.012 (^^^ X 10"^* Wm-^ (2) 
and the line luminosities with: 

^"-^•^^><^°^(i(F^)(fe)^^«' 

where Dz, is the luminosity distance. 
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4.1 SDP.81 

The most prominent features in the FTS spectrum of SDP.81 
are the [Olll] 88 /im and [Cll] 158 /im lines, clearly seen in 
Figure [2l left panel. The continuum subtracted spectral seg- 
ments of SDP.81, at ±5000 kms"^ rest-frame around both 
lines, are shown in Figure O The derived emission line pa- 
rameters from the Sine model fit are given in Table (2] 

We detect the ^P3/2 — ^ P1/2 [Cll] line at 158 ^im, at 
a signal-to-noise ratio (SNR) of ~ 5 (see Figure [21 left 
panel). The line falls near the longer wavelength end of the 
FTS spectrum, in a region in which the instrument con- 
tribution dominates even above the telescope background 
emission and consequently it is visibly noisier (see Figure [21 
left panel). The other prominent feature is the ^Pi — ^ Po 
[O III] 88/im l ine, which is the thir d detection of this line 
at 2 > 0.05 (jFerkinhoff et al.ll201ol l. We do not detect at 
a sufficient significance any of the other possible ISM lines 
that fall in the FTS range. For those lines we give the upper 
limits in Table [3l 



4.2 SDP.130 

We do not have any positive line detection in the FTS spec- 
trum of SDP.130. Even stacking the spectra around all the 
lines in Tables [2] and [3l do not result in any positive signal. 
This target is about 30% fainter than SDP.81 at 350 /im (see 
Table [1] and Figure [2]). Scaling from SDP.81 and assuming 
the same luminosity of the [Cll] 158 /^m line, we estimate 
the expected line flux for SDP.130 to be ~ 4.5 x 10"^** W 
m~^, taking into account the magnification factor. This is 
~ 3 times fainter than the measured integrated line fiux of 
SDP.81. Consequently, to detect this line at the same sig- 
nificance would require 9 times longer integration time. On 
the other hand, the CO (1-0) line lu minosity of SDP.13 is 
50% higher than that of SDP.81 fsee lFraver et al.ll201ll '). so 
scaling directly with the molecular gas mass, we would ex- 
pect to see a marginal, ~ 3-3. 5cr, [C 11] 158 /^m feature. 
There could be diffe rent reasons for a suppressed [C 11] (e.g. 
iLuhman et al .I2OO3I ) but also the scaling of the [C 11] with the 
molecular mass may not be straightforward and tentatively 
may indicate different conditions of the ISM in SDP.130. 

Because of the non-detection we have used the observa- 
tion of SDP.130 as a second dark sky measurement, which 
helped to minimise further the noise in the SDP.81 spectrum 
(see Section [2TT|| . 



4.3 Physical characteristics of the emitting region 
in SDP.81 



We fo llow the analysis p r esent ed in lHailev-Dunsheath et al] 
(|2010l ) and lStacev eTal] (I2010I) using the photo-d issociation 
regions (PDR) models of [Kaufman et al.l (| 19991 ) with the 
help of the PDR toolkit Q The expected physical conditions 
from the observed lines and the total infrared luminosity are 
given in terms of PDR cloud density n in cm^'' and the FUV 
{Q < hv < 13.6 eV) ionization field strength Go, in units of 
the Habing field - the local Galactic interstellar radiation 



field of 1.6 X 10~^ erg cm~^ s~^. We must emphasise that the 
results we derive are for a particular set of PDR models with 
pre-defined parameters, like the visual extinction of Ay = 10 
and metallicity equal to the local ISM metallicity, i.e. Z — 
Zq. At much lower metallicities, i.e. Z = 0.1Z(7) and smaller 







Ay we wo uld expect that the de rived n and Go will change 
by ~ 10% ([Kaufman et al.|[l999l ). In addition, as SDP.81 is 
unresolved, the derived PDR characteristics are an ensemble 
average of possibly many photo-dissociation regions falling 
in the SPIRE-FTS beam. This has an implication f or the 
comparison with the PDR models of [Kaufman et al] (|l999l ) 
using optically thick lines, like the CO transitions or the 
[Ol] 63^tm line. The CO (1-0) emission in SMGs typically 
comes from a relat ively smooth, optically thick medium (e.g. 
[CariUi et aL|[201GD with some substructure that may become 
more prominent in higher-J lines and the continuum. For 
such a medium, CO observations could probe higher by a 
factor of 2 column densities. Nonetheless, in the subsequent 



anal ysis, following [Kaufman et al.l ([l999j) and o ther works 
(e.g. lSwinbank et al.[[2010r Danielson et al.|[2011^ . we double 



the observed luminosities of the optically thick lines, in order 
to account for the fact that we see the emission only from 
one side of the clouds. 



4.3.1 [CllJ 158 fim line 

The [C 11] line at 158 /im is one of the best studied ISM col- 
lisionally excited lines. It is a widespread coolant for the 
ISM in galaxies; it requires only relatively modest FUV 
fiuxes to be produced, and hence it traces star formation 
rather than AGN activity. This line can account for up to 
1% of the total fa r-IR luminosity in normal galaxies (e.g. 
[Stacev et al.[ [l99l[ ) and is optically thin in most cases. A 
large fraction of the [C 11] luminosity comes from photo- 
dissociation regions, however a non-negligible fraction of 
the [C 11] emission may arise from diffuse ionised gas and 
from H II regions. The usual approach to es timate this frac- 
tion is with the [N 11] 122 u m line (see e . g. Malhotra et al 



200 ll: [Contursi et al] [200I [Oberst et all [20061 : [Vasta et al 



2010t ). Nitrogen has ionization potential 14.53 eV, higher 



■* The PDR modelling toolkit llPound fc Wolfirel2007h is available 
at [http : //dustem. astro .umd. edu[ 



than the 13.6 eV of Hydrogen, and so the [N 11] 122 
line can only originate from Hll regions or from the dif- 
fuse ionized medium. When there is not enough informa- 
tion from other lines, which is the case for most of the 
high redshift submm galaxies, it is assumed that 70% of the 
JC III emission arises in PDR (e. g. [Hailev-Dunsheath et al.[ 
[2010l : [Stacev et al.[[201ol ) based on [Oberst et alj (|2006h mea- 
surements in our Galaxy (the Great Carina nebula) and 
also supported by ob servations of nearby normal galaxies 
([Malhotra et al.|[200l[ ). We follow the same approach and 
adopt the same 70% fraction of [C 11] coming from a PDR, 
which is consistent with the observed lower limit of the 
[Cll]/[Nll] ratio of 2.6 in SDP.81, which, compared to the 
ratio of 1.6 from a pure HII region, means that at least 40% 
of the [C 11] flux comes from PDRs. 

The ionising UV source can be either hot stars or an 
AGN and this increases the complexity when one wants 
to derive the physical conditions based on cooling lines in 
the FIR. A line ratio, involving [C 11] , that can be used 
to infer the orig i n of the radiation field, is Lfoii 63^ /-^' [cii] 
([Abel et al.[[2009l : [Meiierink. Spaans fc IsraefcOOT ). The up- 
per limit of this ratio for SDP.81 is 1.1, which is charac- 
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Rest-frame velocity (km/s) Rest-frame velocity (km/s) 

Figure 3. Continuum subtracted regions witiiin ±5000 km s~^ rest-frame for [Cll] 157.7 /jm (left) and [Olll] 88.4/im (right) lines of 

SDP.81. The FTS data processed using two reference dark sky observations (see Section [Sj is shown as histogram. The dashed line 
indicates the zero level. The best fit Sine function is also shown. 



Table 2. Characteristics of the [Cll] and [Olll] lines in SDP.81, as derived with a Sine function fit. We use /iL^lRCS — 2000^m) = 
(5.13 ± 0.88) X 10^ ^ LfTi and = 25± 7, based on an improved lens model for SDP.81, which is different but still compatible with 
/^L = 19 quoted in lNegrello et all l l2010t) . The reported velocity offset vq, in rest-frame kms~^, is between the fitted line centroid and 
the expected one for the CO measured redshift z = 3.042. The instrumental line FWHM in kms^^ corresponds to the high spectral 
resolution of 1.2 GHz, we keep this parameter fixed for the Sine fit. The noise rms is the standard deviation of the residual (fitted model 
subtracted from the data) in the velocity range ±5000 kms~^. The reported signal-to-noise ratio (SNR) is the fitted line peak divided 
by the noise rms. The integrated line flux and line luminosity are calculated using Eqs. [2] and |3] respectively. The errors are statistical 
errors, they do not include the FTS flux calibration uncertainty, conservatively estimated to be 15-20% for wavelengths less than 300 
lira and 20-30% for wavelengths above 300 ^m. 



Characteristics [Olll] 88.4 /im [Cll] 157.7 /^m 

Fitted line centroid Ac {fira) 357.2 ± 0.6 (839 GHz) 638.4 ± 1.8 (490 GHz) 



Redshift 


3.041 ± 0.006 


3.046 ±0.011 


Velocity offset iiq (kms~^) 


95 ±30 


-240 ± 60 


Instrumental line FWHM fpwHM (kms"'^) 


514 


919 


Fitted line peak /(O) (mjy) 


480 ± 70 


1100 ± 160 


Noise rms (mJy) 


90 


211 


SNR 


5.4 


5.2 


Line flux Sune (xlO"!* W m'^) 


5.8 ±0.8 


13.2 ± 1.9 


Velocity integrated line flux Sn^c (Jy kms^^) 


210 ±30 


840 ± 120 


Line luminosity /iL^iinc (xlO^" Lq) 


13 ±2 


29 ±4 


ilino/ilR (xlO-3) 


2.4 ±0.5 


5.6 ±1.3 



teristic for environments excited by stellar radiation, but 
also overlaps marginally with the allowed region for AGN 
excitation, which would require a line ratio greater than 1 
I Abel et al.ll2009l ). Thus, we have to keep in mind that part 
of the ionizing radiation field may have contribution from an 
AGN. T he FIR lines are not sensit ive indicator of an AGN 
(see e.g. iGracia-Carpio et al.ll201ll ) which makes it difficult 
to quantify this fraction. In the subsequent PDR analysis we 
assume that the radiation field is due to UV coming from 
hot stars. 

An idea of the physical conditions, in terms of the 
FUV radiation field strength Go and PDR density n, can 
be obtained by comparing the fractional [C ll] luminos- 
ity vers us the fractional CO(l-O) luminosity with PDR 
models (jStacev et all Il99ll: iHailey-Dunsheath et al.1 l2010l : 
llvison et al. I l2010dl : IStacev et al.l |2010| ). In the diagram 
shown in Figure [4l where we have used the EVLA C0(1- 



Table 3. Upper limits (3(t) for some ISM lines for SDP.81 that 
fall in the FTS spectral range. The upper limits were calculated 
using Eq. |2] assuming /(O) to be 3(7 of the noise st.dev. within 
±5000 kms~'^ around the expected line centre. 



Line ID 


Ao 


■^linc 








(^tm) 


(10-18 W m-2) 


(lOi°i0) 


(xlO-3) 


[Oiii] 


52.0 


< 4 


< 9 


< 1.7 


[N III] 


57.0 


< 4 


< 10 


< 1.9 


[Oi] 


63.2 


< 5 


< 11 


< 2.2 


[Nil] 


122.1 


< 5 


< 11 


< 2.2 


[Oi] 


145.5 


< 8 


< 19 


< 3.6 



0) measurements from this work, we can see that even with 
the [C ll] alone SDP.81 is comprised of PDRs with a mean 
density of ~ 10^ - lO" cm"^ and Go ~ 10^ - 10^, similar to 
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other star-forminR g alaxies at 2; = 1 — 2 and to SMMJ2135 
lllvison et al ] |2010dh . We can compare the relative strength 
of the [C 11] line, in terms of the FIR lumino sity, with other 
galaxies with measured [Cll] lines (see e.g 



Maiolino et alj 



201 ll ). SDP.81 



I2OO9I : IStacev et aLlbOld : iGracia-Carpio et al 
has log (L[cji]/Lfir) = —2.2 ± 0.05 - one of the highest ob- 
served ratios, which contra sts with the deficien cy reported 
for ULIRG and AGN (e.g. iLuhman et al.llioo^ '). Although 
it might seem that SDP.81 relative [C 11] strength is much 
higher than the one seen in AGNs at 2 = 1 — 2, recent mea- 
surements with the Herschel FAGS instrument, however, do 
not show any clear separation between local AGN systems 
(QSO, S eyferts, LINERS), starburs ts and infrared luminous 
galaxies (|Gracia-Garpio et al.ll201ll ). 

The ratio i[cii]/^co(i-o) of ~ 4000 is quite close to the 
average of 4100 for normal metallicity star-formi ng galax- 
ies or Galactic star-form ation or molecular clouds (|Maddenl 
l2000l : IStacev et al.ll2010D . It is also important to note that in 
comparison with the other high redshift galaxies SDP.81 has 
one of the highest fractional [Gll] and CO(l-O) luminosities. 
Whether this is significant or it is simply due to luminosity 
bias is to be confirmed by increasing the number of high 
redshift SMGs at Lfir < 10^^ Lq. 

The diagram shown in Figures U is quite useful as diag- 
nostic and as a comparison with other observations of galax- 
ies of different classes. To properly quantify the physical con- 
ditions, we can add all the available information: the FDR 
luminosity of [G 11] , the total FIR luminosity, the upper limit 
of the [Ol] 63/im hue and the EVLA detections of GO (1-0) 
and perform a least square fit with grids of FDR models with 
differ ent (n. Go). Using the FDR toolkit (|Found fc Wolfird 
I2OO7I ) we derive the best fit FDR density of n = io3-3''±o.22 
cm~'^, which is in very good agreement wi th the H2 density 
from RADEX models using the GO lines (|Lupu et al.ll2010l . 
model m03). The corresponding radiation field strength is 
Go = io2-2-4±o.i9^ rpj^g confidence region in the (n,Go) 
plane, delineated by the line ratios, is shown in Figure [S] 
These (n, Go) values correspond to a FDR surface temper- 
ature of Ts ~ 200 K. 

Adding the GO(7-6) hue from ILupu et all (|2010l ). and 
considering only the two ratios [G ll]/GO(l-0) and GO(7- 
6) /CO (1-0), we obtain an order of magnitude higher val- 
ues: n ~ 3.2 X 10* cm"'' and Go ~ 5.6 x 10^, as shown in 
Figure [5] This is to be expected as the GO(7-6) line usu- 
ally traces higher density regions. In addition, it must be 
noted that GO(7-6) is likely to be blended with [Cl] given 
the relatively coarse resolution of the Z-Spec spectrometer. 
Recent measurements with Herschel show that [C l] has ap- 
proximately 50-150% of the flux of the GO(7-6) line, 80% 
on average, based on a sample of 20 local ultra-luminous 
infrared galaxies (van der Werf et al, in preparation). Even 
taking into account the blending, i.e. scaling down the G0(7- 
6) line flux by a factor of 1.8, the overall picture does not 
change significantly. It should be noted that Go derived from 
[G 11] should be more representative for the FDR as a whole, 
since the [C 11] line comes from FDRs and molecular cloud 
surfaces. Gonsequently, in the following analysis we use the 
low (n. Go) values (i.e. not taking into account the GO(7-6) 
line). 

A straightforward use of the FDR derived characteris- 
tics is to investigate the effective physical size of the non- 
lensed emitting region. With FDR density n of ~ 2000 



cm and the total gas mass from ILupu et all (|2010l ) of 
Mgas = 1.3 X IO^^A/qZ/il, where fii, is the magnification 
factor of 25 ± 7, we obtain the effective radius of the FDR 
emitting region of ~ 1.0 kpc/fi]/^ or (0.3-0.4) kpc for a 
spherical distribution. For a disk geometry with thickness 
of 100 pc, the radius is ~ 3.0 kpc/y^T'L or (0.5-0.7) kpc 
correspondingly. This result is in good agreement with the 
eq uivalent dus t emitt ing radius of 0.7 kpc for SDF.81 derived 
in lLupu et a .1 (|2010|). and the size of the emitting region in 
SMMJ2135 (|Danielson et al.ll201ll ). 

Although the SFIRE FTS smallest beam size is 17" 
FWHM, nevertheless we can have an idea of the the effec- 
tive magnified size of the [G 11] emitting region, by means of 
the [G 11] beam filling factor. Assuming that the FDR slab 
is illuminated from one side by a flux Go and that dust 
grains are heated by the FUV photons th en we can com- 
pare the model predicted [G 11] line intensity (|Kaufman et al.l 
I1999I ) with the observe d line intensity and thus we can derive 
the beam fiUing factor (IWolfire. Tielens fc Hollenbach|[l990l ) 
*A = (^[ciii/fi) /I[cu] « 3 X 10"*, where n = 10"'' sr = 
1.22 arcmin^ is the SFIRE FTS effective b eam area at the 
wavelength of the [G 11] line (|Swinvard et al.i r201Q). This cor- 
responds to an effective radius of the magnified emitting re- 
gion of « 1". This is consistent with the SMA continuum 
observations and the EVLA GO(l-O) map, shown in Fig- 
ure [T] 

We can also use the FDR results to place a lower 
limit on the H2 mass of SDF.81, independently from 
the mass derived using GO measurements. Assuming that 
the [G 11] line is optically thin the n the minimum total 
atomic hydrogen mass in the FDR (iGrawford et all Il985l : 
iHailev-Dunsheath et al.ll2010l : iKdson et al.ll2010dl ): 



M. = O.77Q(r)^Hiii^M0, 



(4) 



^0 X(C+) 

where Q{T) = [l + 2e("«i + n,„,/n] / [2e^-''' ^/^-')] . 

Substituting the [G 11] luminosity corrected for the H 11 con- 
tribution and for the magnification, the FDR surface tem- 
perature and density, the critical density of G "*" '^P2/3 level 
of Ucrit = 2700 cm"^ (iLaunav fc Rouefall977l) a nd the G+ 
abundance x(C+) = 1.4 x 10 * (jSavage fc SembachT996i ). 
we derive Ma « 1.7 x 10^°Mq. It should be noted that this 
minimum atomic mass in the FDR is quite sensitive to the 
assumed FDR fraction of the [G 11] luminosity, the magni- 
fication factor and the unknown real heavy element abun- 
dances. Nevertheless the derived Ma is of the order of the 
estimated molecular gas mass M(H2) = (1 - 3) X 10^"^ Mq 
derived from GO(l-O) (jPraver et al.ll201 j ), somewhat higher 
than the fraction of 20-50% seen in o ther high-redshift star- 
forming galaxies jStacev et al.|[201oh . 

4.3.2 [Oni] 88 ^im line 

The [O III] lines are usually connected with H 11 regions but 
also with AGN activity. The [O ill] 88fim line is emitted via 
the transition ^Pi — ^ Po of O^""". The ionization potential 
of O^ is 35 eV, which requires a nearby hot young star or 
stars (T ~ 36 000 - 40 000 K) or an AGN, to provide the 
necessary FUV photons. The [O ill] line emission is quite 
sensitive to the effective temperature of the ionizing stars 
and if this temperature drops from 36,000 K to 33,000 K 



10 Valtchanov et al. 




Figure 4. Diagnostic diagram of L[qji]/^fir as a function of 1/00(1— 0)/^FlR f°'' SDP.81 (the rightmost fi lled circle) an d SDP.130 
(filled circle with a down arro w) , compared to local and other known high redshift galaxies of different c lasses j Stacey et al.||l99ll. 120101 : 
iHailev-Dunsheath et aLll2010l ). The lines with constant density n and constant Gq from PDR models of iKaufinaii^taLrilQQsl ) are also 
shown. The [Cll] luminosity was corrected for the fraction of the flux arising from Hll and diffuse ionized regions. The CO(l-O) line 
flux from the EVLA measurement was multiplied by 2 to take into account that the line is optically thick and we see only the fraction 
coming from the illuminated PDR side. i[cil]/^CO(l-0) of 4100, typical for local star-burst galaxies, is also shown as a dashed line. We 
have used Lfir(42.5 - 122.5Aim) = Lir(3 - 2000/xm)/1.54, derived by integrating the SED of SDP.81. 



the li ne intensity would drop by a factor of up to 30 jRubinl 
1 19851 ). The ratio L[oiii]/-^'Fir observed in SDP.81 is signifi- 
cantly higher than the one observed in normal galaxies, al- 
though s maller than the ratio in irregular, Hll dominate d 
galaxies (|Malhotra et al.ll200ll : iGracia-Carpio et aLll2011^ . 
Comparing this ratio with models ( Abel et al.ll2009l ) indicate 
that some fraction of the emission may be due to an AGN 
source. The situation is not that clear, however, comparing 
L[oiii]/ Lfik with recent observations, which show no clear 
separation of the l ocal AGN sources (Seyfer ts, QSOs, LIN- 
ERs) and the rest (|Gracia-Carpio et al.ll20l"ll ). Thus, similar 
to the ^[oi] 63/j/^[c II] result, we cannot discard completely 
that some fraction of the ionizing field may arise in an AGN. 

The detected [O ill] 88 ^ra line can also provide us with 
information about the properties of an ensemble H ll compo- 
nent, with [Oiii] emission coming from region s much closer 
to th e UV sources (see e.g. I Cormier et al ] |2010l : l0kada et al.l 
l20ld ). Transitions of the same ion O"*"^ from different col- 
lisionally excited levels of the same fine-structure multiplet, 
as is the case for the two lines of [Om] at 52 and 88 /im, 
are not sensitive to the temperature, for a reasonably large 
range of the density, and can be used to estimate the elec- 



tron density rie in the H ll region. Unfortunately we do not 
detect [O ill] at 52 /im, thus we can only put an upper limit 
on i[oiii]52M/^[oiii]88M of ~ 0-7, which is quite c lose to 
the low density limit of Tie ~ 100 cm~'^ (see e.g. iRubinI 
1 19891 : iFischer et al] Il999l ). This electron density is typical 
for galactic Hll regions. 



4.4 FIR/radio correlation and millimeter imaging 
source constraints 

In support for the hypothesis that some fraction of the ion- 
izing radiation in SDP.81 may originate from an AGN is 
the radio data (see Section l2.3p . In Figure [6] we compare 
the Spe ctral Energy Distribu tion of SDP.81 to the template 
used in llvison et all (|2010dl ) for the SMMJ2135 (the cos- 
mic "Eyelash"), which follows the FIR/radio correlation. 
Combining 5'8. 4GHz with 5'i.4GHz = 950 ± 150 /iJy from the 
FIRST survey (|Becker et al.|[l9 95l) yields a spectral index of 
Q = — 0.17to'o9 (where oc v"), suggesting th at an AGN 
is pow ering some fraction of the radio emission (|lbar et al.l 
l2010al ). Note however that a flat spectrum could also be 
a result from compact star-forming regions dominated by 
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Figure 5. PDR model predictions for the density (n) and the 
FUV fiux (Go) from the observed line ratios. The [O l] 63/im/[C ll] 
158 /im is shown as an upper limit, the dashed line indicating the 
side where we expect the ratio to be if [O l] is less luminous. The 
shaded contours show the 3(t confidence region of the expected 
(n, Go) values when considering only [C ll], CO(l-O) and the upper 
hmit on [Ol]. The CO(7-6)/CO(l-0) ratio was corrected for [Cl] 
blending with CO{7-6). 



star formation l|Hunt fc Maiolin"oll2005^ ■ Combining the ra- 
dio measurements with Lir yields gia = 1.8 ± 0.1 (see 
llvison et al.ll20103 for details) . This is only barely consistent 
with the FIR-radio correlation, providing more evidence 
that SDP.81 is moderat ely radio-loud (cf. g iR = 2.4 for FIR- 
selected galaxies - see [ivison et al.l I2010al l31dl: Ijarvis et al.l 
l2010l : lMichalowski et al.ll2010l ). 

At mm-band, the peak and the total 1 mm fluxes mea- 
sured in a PdBI 0.62" x 0.30" beam and ~ 3" diameter 
area of SDP.81 are 2.4 mjy and 7.5 mjy respectively (Neri 
et al. in preparation), suggesting that the contribution of 
any strong compact (< 0.54") source, like an AGN - if it is 
indeed present in this object, would be at most 33% of the 
total mm radiation. If the unresolved 20 mJy flux detected 
by the TRAM 30 m telescope at 2 mm is also considered 
jNegrello et al.ll2010l ). then the PdBI measured peak sug- 
gests the contribution of any strong compact source would 
be less than 14% of the total flux at these wavelengths. 

With higher resolution observations using CARMA in 
its most extended configuration, SDP.81 is not detected. The 
rms noise in the 0.17" x 0.14" beam map is estimated at 
0.67 mJy. Based on this and the peak emission in the PdBI 
observations, we conclude that if the peak emission in the 
PdBI image was a point source, it would have been detected; 
therefore, it has to be at least 0.15" to be consistent with 
the CARMA non-detection. This impUes the 2.4 mJy peak 
emission measured in the PdBI beam cannot all come from 
a point source or AGN, putting an upper limit of the AGN 
flux contribution to no more than ~ 1/3 (or 0.8 mJy) of the 
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Figure 6. SED of SDP.81 including the new measurements pre- 
sented in this paper: the radio observation at 8.5 GHz (see Sec- 
tion l2.3l l and PACS 70 and 160 fim, all marked as di amonds. The 
1.4 G Hz detection of SDP.81 in the FIRST survey l|Becker et al.l 
Il995h is also shown. The dashed line is an Arp 220 template, as 
used for SMMJ2135 source (the cosmic "Eyelash") reported in 
llvison et al.l l l2010dl ) and that obeys the FIR/radio correla tion. 
The filled circles repr esent data from iNegrello et al ] l|2010h and 
iHopwood et al.l 1 I2OIII ). The two radio measurements are clearly 
deviating from the FIR/radio correlation followed by other FIR 
selected galaxies. 



2.4 mJy peak and ~ 1/10 of the 7.5 mJy total mm emis- 
sion detected by PdBI or ~ 1/25 of the 20mJy detected by 
the IRAM telescope. This is consistent with the continuum 
upper limit of 1 mJy measured at 1 cm in the EVLA obser- 
vations, when extrapolated to the 1 mm wavelength regime 
for a radio-loud or flat-spectrum source. 

Putting all the pieces together we can see that the 
physical conditions in SDP.81 are similar to other high red- 
shift star-forming galaxies, although SDP.81 is at somewhat 
lower intrinsic luminosity, but still being an ultra-luminous 
infrared galaxy. The radio observations provide supporting 
evidence to the hints from the ISM lines, that part the ion- 
ization field in the galaxy may be due to an AGN. With the 
available observations we cannot put any strong constraint 
on this fraction, although a compact source contribution in 
the mm-band fiux cannot exceed 33%. 



5 CONCLUSIONS 

Wc have presented //er5c/ie/-SPIRE Fourier Transform 
Spectrometer and radio follow-up observations of Herschel- 
ATLAS strongly lensed sources. Two targets, SDP.81 and 
SDP.130 were observed with the FTS, resulting in spectra 
over a wavelength range from 197 to 670 pm (447-1522 GHz) 
at a constant spectral resolution of 1.2 GHz. We also pre- 
sented observations of SDP.81 with the EVLA, PdBI and 
CARMA in the radio domain, helping to better characterise 
the source properties. 

Thanks to the lensing effect, boosting the luminosity of 
SDP.81 by a factor of ~ 25, we were able to detect [Cll] at 
158 fim and [O ill] 88 ^m lines in its FTS spectra. We do not 
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detect any plausible lines in the spectrum of SDP.130. We 
should have been able to marginally detect the [C ii] 158 ^m 
line at 3-3. Sct if we scale the [C ll] with the molecular mass 
derived from ground-based CO(l-O) line measurements. If 
we assume the same ratio of [Cll] to Lm as in SDP.81, how- 
ever, we do not expect detection - the line would be too 
faint. This is an indication of a possible difference in the 
quantity and the distribution of the molecular material in 
the two FTS targets, although the uncertainties are too im- 
portant for any strong claim. 

The successful line detections in SDP.81 allowed for de- 
tailed studies of the physical conditions in the emitting re- 
gions at 2 ~ 3. The clear and unambiguous detections of 
[O III] 88^m and [C ll] 158 /^m allow to measure the redshift 
using both lines. The derived value of 2 = 3.043 ±0.012, us- 
ing the two lines, confirms the measured redshift of 2 = 3.042 
from ground based CO measurements. The detection of the 
[O III] 88 ^m line is the third at redshifts higher than 0.05. 
We expect that the on-going SPIRE and PACS spectroscopy 
with Herschel will lead to many more such line detections. 

Combining the FTS detected [C 11] 158 nm line, the total 
infrared luminosity, the upper limit on the [O l] 63 fim line, 
the ground based detections of the CO(l-O) transition we 
were able to infer a PDR cloud-ensemble density no ~ 2000 
cm~"^ and the far-UV ionizing field strength Go ~ 200 in 
SDP.81. The PDR surface temperature is about 200 K and 
the effective radius of the PDR emitting region is of the 
order of 500-700 pc for a disk of 100 pc thickness. These 
characteristics are similar to other high redshift star-forming 
galaxies. 

The presence of a strong [O ill] line, the upper limit 
of i[oi]63/j/^[cii] and the significant deviation from the 
FIR/radio correlation are all plausible indications that some 
fraction of the emission from SDP.81 may be due to an AGN 
source. The conservative upper limit of this contribution in 
the radio domain is estimated at ~ 33 %. 

The analysis presented in this paper is a clear illustra- 
tion of the necessity of a combined effort from follow-up ob- 
servations at different wavelengths in order to have a better 
idea of the physical conditions in the star- forming galaxies in 
the distant Universe. Herschel will be indispensable in this 
task, providing large samples of such galaxies for follow-up 
studies. 
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